Kinematic entanglement degradation of fermionic cavity modes 
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We analyse the entanglement and the nonlocality of a (1 + l)-dimensional massless Dirac field 
confined to a cavity on a worldtube that consists of inertial and uniformly accelerated segments, for 
small accelerations but arbitrarily long travel times. The correlations between the accelerated field 
modes and the modes in an inertial reference cavity are periodic in the durations of the individual 
trajectory segments, and degradation of the correlations can be entirely avoided by fine-tuning 
the individual or relative durations of the segments. Analytic results for selected trajectories are 
presented. Differences from the corresponding bosonic correlations are identified and extensions to 
massive fermions are discussed. 
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I. INTRODUCTION 

One of the fundamental problems in the emerging field 
of relativistic quantum information is the degradation of 
correlations caused by accelerated motion. Studies of uni- 
form acceleration in Minkowski spacetime (see p]-[(| for 
a small selection and Q for a recent review) have re- 
vealed significant differences in the degradation that oc- 
curs for bosonic and fermionic fields. There are in par- 
ticular clear qualitative differences in the bosonic versus 
fermionic particle-antiparticle entanglement swapping Q 
and in the infinite acceleration residual entanglement and 
nonlocality 

The analyses of uniform acceleration mentioned above 
involve two ingredients that make it difficult to com- 
pare the theoretical predictions to experimentally real- 
isable situations. The first is that while the uniformly- 
accelerated observers are considered to be pointlike and 
perfectly localised on a trajectory of prescribed accelera- 
tion, the field excitations are nevertheless usually treated 
as delocalised field modes of plane wave type, normalised 
in the sense of Dirac rather than Kronecker deltas. This 
may seem a technicality-perhaps remediable by use of ap- 
propriate wave packets [J| , but at present it appears unex- 
plored how localised observers would in practice perform 
measurements to probe the correlations in the delocalised 
states. 

The second concern lies in the time evolution of the 
correlations. An inertial trajectory in Minkowski space is 
stationary, in the sense that it is the integral curve of a 
Minkowski time translation Killing vector. A uniformly- 
accelerated trajectory is also stationary, in the sense that 
it is the integral curve of a boost Killing vector. How- 
ever, the combined system of the two trajectories is not 
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stationary, as the two Killing vectors do not commute. 
For example, in the (1 + l)-dimensional setting there is a 
unique moment at which the two trajectories are parallel, 
and the trajectories may or may not intersect depending 
on their relative spatial location. Yet the analyses men- 
tioned above regard the correlations between observers 
on the two trajectories as stationary and the relative lo- 
cation of the trajectories as irrelevant, observing just that 
the spacetime has a quadrant causally disconnected from 
the uniformly-accelerated worldline and noting that the 
field modes confined in this quadrant are inaccessible to 
the accelerated observer. While the acceleration horizon 
that is responsible for this inaccessibility may be seen as 
the basis of the Unruh effect [8, 9] , the horizon exists only 
if the uniform acceleration persists from the asymptotic 
past to the asymptotic future. In this setting it is not 
clear how to address motion on trajectories that remain 
uniformly accelerated only up to the moment at which lo- 
calised observers might make their measurements on the 
quantum state. 

Both of these concerns have been recently addressed by 
studying correlations between field modes of a real scalar 
field confined in two cavities, one inertial and the other 
undergoing motion that consists of segments of inertial 
motion and uniform acceleration (Toj . In this setting the 
field modes are spatially localised in the cavities, and 
the acceleration effects can be localised in time by taking 
the initial and final segments of the accelerated cavity 
to be inertial. It was found that the entanglement is 
affected by the acceleration, and in (1 + l)-dimensional 
spacetime the mass of the field has a strong effect on 
the qualitative behaviour on the entanglement. For a 
massless Dirichlet field the entanglement is periodic in 
the durations of the individual trajectory segments, so 
that entanglement degradation can be entirely avoided 
by fine-tuning the durations of the individual segments; 
further, in the small acceleration limit the degradation 
can also be avoided by fine-tuning the relative lengths of 
the inertial and accelerated segments. 

In this paper we shall undertake the first steps of in- 
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vestigating fermionic entanglement in accelerated cavities 
by adapting the scalar field analysis of [h| to a Dirac 
fermion. Conceptually, one new issue with fermions is 
that the presence of positive and negative charges allows 
a broader range of initial Bell-type states to be consid- 
ered. Another conceptual issue is that in a fermionic Fock 
space the entanglement between the cavities can be char- 
acterised not just by the negativity but also by the viola- 
tion of the Clauser-Horne-Shimony-Holt (CHSH) version 
of Bell's inequality [ll], [HI, physically interpretable as 
nonlocality. New technical issues arise from the bound- 
ary conditions that are required to keep the fermionic 
field confined in the cavities. 

We focus this paper on a massless fermion in (1 + f) 
dimensions. In this setting another new technical issue 
arises from a zero mode that is present in the cavity under 
boundary conditions that may be considered physically 
preferred. This zero mode needs to be regularised in order 
to apply usual Fock space techniques. 

We shall find that the entanglement behaviour of the 
massless Dirac fermion is broadly similar to that found 
for the massless scalar in [Toj , in particular in the periodic 
dependence of the entanglement on the durations of the 
individual accelerated and inertial segments, and in the 
property that entanglement degradation caused by accel- 
erated segments can be cancelled in the leading order in 
the small acceleration expansion by fine-tuning the dura- 
tions of the inertial segments. We shall however find that 
the charge of the fermionic excitations has a quantita- 
tive effect on the entanglement, and there is in particular 
interference between excitations of opposite charge. 

We begin in Sec.|II]by quantising a massless Dirac field 
in a static cavity and in a uniformly-accelerated cavity in 
(I + 1) dimensional Minkowski spacetime. We pay special 
attention to the boundary conditions that are required for 
maintaining unitarity and to the regularisation of a zero 
mode that arises under an arguably natural choice of the 
boundary conditions. Section IIIII develops the Bogoli- 
ubov transformation technique for grafting inertial and 
uniformly accelerated trajectory segments, presenting the 
general building block formalism and giving detailed re- 
sults for a trajectory where initial and final inertial seg- 
ments are joined by one uniformly accelerated segment. 
The evolution of initially maximally entangled states is 
analysed in Sec. IIV1 and the results for entanglement 
are presented in Sec. [V] A one- way-trip travel scenario, 
in which the accelerated cavity undergoes both accelera- 
tion and deceleration, is analysed in SecEU Section IVTll 
presents a brief discussion and concluding remarks. 

We use units in which h = c = 1. Complex conjugation 
is denoted by an asterisk and Hermitian conjugation by 
a dagger. 0(x) denotes a quantity for which 0{x)/x is 
bounded as x — > 0. 



II. STATIC CAVITY 

In this section we quantise the massless Dirac field in 
an inertial cavity and in a uniformly accelerated cavity, 
establishing the notation and conventions for use in the 
later sections. 



A. Inertial cavity 

Let (£, z) be standard Minkowski coordinates in (1 + I) 
dimensional Minkowski space, and let rj^ denote the 
Minkowski metric, ds 2 = dx^ dx u = —dt 2 + dz 2 . The 
massless Dirac equation reads 

i^d^=0, (1) 

where the 4x4 matrices 7^ form the usual Dirac-Clifford 
algebra, {7^,7"} = 2j] fJ - l/ . A standard basis of plane 
wave solutions reads 

^W(t, z) = A w , e , CT e"^'-^ U e , a , (2) 

where u> G R, e G {1,— I}, a G {1,-1}, the constant 
spinors U e>(T satisfy 

OfiUe,a = et^cr , (3a) 

7 5 f/ e , (T = aU c . a , (3b) 
UlrUw = <W<W< , (3c) 

and A^.g.o- is a normalisation constant. Physically, uj 
is the frequency with respect to the Minkowski time, 
the eigenvalue e of the operator = 7°7 3 indicates 
whether the solution is a right-mover (e = 1) or a left- 
mover (e = —1), and a is the eigenvalue of the helic- 
ity/chirality operator 7 s = i7°7 1 7 2 7 3 [l^]. The right- 
handed (a = +1) and left-handed (a = —1) solutions are 
decoupled because flU does not contain a mass term. 

We encase the field in the inertial cavity a < z < b, 
where a and b arc positive parameters satisfying a < b. 
The inner product reads 

6 

(i>(i),i>(2)) = J dzip\ 1) ip( 2 ) , (4) 

a 

where the integral is evaluated on a surface of constant t. 
To ensure unitarity of the time evolution, so that the 
inner product Q is conserved in time, the Hamiltonian 
must be defined as a self-adjoint operator by introducing 
suitable boundary conditions at z = a and z = b fLU [l5| . 
We specialise to boundary conditions that do not couple 
right-handed and left-handed spinors. For concreteness, 
we consider from now on only left-handed spinors, and we 
drop the index a. The analysis for right-handed spinors 
is similar. 
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We seek the eigenfunctions of the Hamiltonian in the 
form 

V> w (i, z) = A u e-^'- 2 ) U+ + e --(*+ z ) JJ_ , (5) 

where A u and B u are complex-valued constants. It would 
be mathematically possible to maintain unitarity by al- 
lowing probability to flow out through one of the cav- 
ity walls and instantaneously reappear at the other wall; 
physically, this would mean that the spatial surface is 
considered to be a circle, possibly with one marked point. 
However, we wish to regard the spatial surface as a gen- 
uine interval with two spatially separated endpoints, and 
we hence specialise to boundary conditions that ensure 
vanishing of the probability current independently at 
each wall. The boundary condition on the eigenfunctions 
thus reads 



(6) 



where ip = ip f-y . 

Following the procedure of [H], [HI , we find from ([5| 
and ([6]) that the self-adjoint extensions of the Hamilto- 
nian are specified by two independent phases, character- 
ising the phase shifts of reflection from the two walls. We 
encode these phases in the parameters <E [0, 27r) and 
s G [0,1), which can be understood respectively as the 
normalised sum and difference of the two phases. The 
quantum theories then fall into two qualitatively differ- 
ent cases, the generic case < s < 1 and the special case 
s = 0. 

In the generic case < s < 1, the orthonormal eigen- 
functions are 

e -iuj„,(t-z+a) jj I £ i0 e -iuj„,(t+z-a) jj _ 

Mt, z) = , (7a) 



■2<) 



(n + s)tt 



(7b) 



where n 6 Z and 8 := b — a. Note that ui n ^ for 
all n, and positive (respectively negative) frequencies are 
obtained for n > {n < 0). A Fock space quantisation 
can be performed in a standard manner [l3j . 

The special case s = corresponds to assuming that 
the two walls are of identical physical build. In this case 
uj n 7^ for n ^ but luq = 0. It follows that a Fock 
quantisation can proceed as usual for the n ^ modes, 
but n = is a zero mode that does not admit a Fock 
space quantisation. In what follows we consider the s = 
quantum theory to be defined by first quantising with 
s > and at the end taking the limit s — > 0+. All our 
entanglement measures will be seen to remain well defined 
in this limit. 



B. Uniformly accelerated cavity 



We consider a cavity whose ends move on the worldlincs 



as above The proper accelerations of the ends arc 1/a 
and 1/b respectively, and the cavity as a whole is static 
in the sense that it is dragged along the boost Killing 
vector £ := zdt + td z . At t = the accelerated cavity 
overlaps precisely with the inertial cavity of Sec. Ill Al 

Coordinates convenient for the accelerated cavity are 
the Rindler coordinates (r],x), defined in the quadrant 
z > \t\ by 



t = x sinh ?/ , z = x cosh rj 



(8) 



where < x < oo and — oo < ?y < oo. The metric reads 
ds 2 = —x^drj 2 + dx 2 ■ The cavity is at a < x < b, and 
the boost Killing vector along which the cavity is dragged 
takes the form £ = d v . 

In Rindler coordinates the massless Dirac equation (TTJ) 
takes the form [l6|, [TtJ 

idv^iViX) = {-ia 3 {xd x + |)) ip(v,X) , (9) 

and the inner product for a field encased in the acceler- 
ated cavity reads 



("0(1)>(2)) = / dx^\ 1) ip(2) 



(10) 



where the integral is evaluated on a surface of constant 77. 
Working as in Sec. Ill Al we find that the orthonormal 
energy eigenfunctions are 



$n(v,x) 



y/2 X HWa) 



(n + s)ir 
ln(6/a) ' 



(11a) 
(lib) 



V a 2 + t 2 and 



\/b 2 + 1 2 , where the notation is 



where n 6 Z. The parameters 9 and s have the same 
meaning and values as above: we consider the micro- 
physical build of the cavity walls not to be affected by 
their acceleration. For s / a Fock space quantisation 
can be performed in a standard manner. For s = the 
mode n = is again a zero mode, and we consider the 
s = quantum theory to be defined as the limit s — > 0+. 



III. GRAFTING TRAJECTORY SEGMENTS 

We now turn to a cavity whose trajectory consists of 
inertial and uniformly accelerated segments. 

The prototype cavity configuration is shown in Fig. [TJ 
Two cavities, referred to as Alice and Rob, are initially 
inertial and in the configuration described in Sec. Ill Al 
At t = 0, Rob's cavity begins to accelerate to the right, 
following the Killing vector £ — d v . The acceleration 
ends at Rindler time r\ — 771, and the duration of the 
acceleration in proper time measured at the centre of the 
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FIG. 1. Space-time diagram of cavity motion is shown. Rob's 
cavity is at rest initially (Region I), then undergoes a period 
of uniform acceleration from t — to r\ — rji (Region H) and 
is thereafter again inertial (Region m). Alice's cavity overlaps 
with Rob's cavity in Region I and remains inertial throughout. 



cavity is ti := | (a + 6)771. We refer to the three segments 
of Rob's trajectory as Regions I, I and IE. Alice remains 
inertial throughout. 

We shall discuss the evolution in Rob's cavity in two 
steps: first from Region I to Region I and then from 
Region I to Region IE. We then use the evolution to 
relate the operators and the vacuum of Region I to those 
in Region H. 



A. Region I to Region II 

Consider the Dirac field in Rob's cavity. In Regions 
I and I we may expand the field using the solutions (|7|l 
and (|TT|) respectively as 



I : tp = a„ i' n + b n tp n , (12a) 

n>0 n<0 

I : ip = ^2 a m i' m + bin , (12b) 

m>0 m<0 

where the nonvanishing anticommutators are 

I : {a m , a n } = {b m , b n } = S mn , (13a) 
I : [a m , a-l} = {6 m , 6* } = 5 mn . (13b) 

Matching the expansions (fT2"|) at t = 0, we have the Bo- 
goliubov transformation 



ipm — ^ ' A mn ip n , ip n — S ' A mn xp n 



(14) 



where the elements of the Bogoliubov coefficient matrix 
A = (Amn) are given by 



A, 



(15) 



and the inner product in (|15[) is evaluated on the surface 
t = 0. Note that A is unitary by construction. 

We shall be working pcrturbatively in the limit where 
the acceleration of Rob's cavity is small. To implement 
this, we follow [l(| and introduce the dimensionless pa- 
rameter h := 25/(a + b), satisfying < h < 2. Physically, 
h is the product of the cavity's length 5 and the accel- 
eration at the centre of the cavity. Expanding (|15[) in a 
Maclaurin scries in h, we find 



A = A<® + + A^ + 0(h 3 ) , 



(16) 



where the superscript indicates the power of h and the 
explicit expressions for A^°\ A^ and A^ read 



-^mn — u mn ; 



= 

1} _ [{-l) m+n - l](m + n + 2s) 



(17a) 
(17b) 



A' 



nn 



2ir 2 (m — n) 3 
1 TT 2 (n + s) 2 



h , (m 7^ n) 



96 



240 



(17c) 
(17d) 



^= tC-n)' ] [(m + S)2 + 3(rt + S)2 

+ 8(m + s)(n + s)]/i 2 . (m ^ n) (17e) 

The expressions (|T7|) show that the small /i expansion of 
A mn is not uniform in the indices, but we have verified 
that the expansion maintains the unitarity of A pcrturba- 
tively to order h 2 and the products of the order h matrices 
in the unitarity identities are unconditionally convergent. 

The perturbative unitarity of A persists in the limit 
s^0|. Had we set s = at the outset and dropped the 
zero mode from the system by hand, the resulting trun- 
cated A would not be perturbatively unitary to order h 2 . 



B. Region I to Region HI 

In Region IH, wc expand the Dirac field in Rob's cavity 



as 



I : V = a » ^« + 12 i>n , (18) 



n>0 



n<0 



where the mode functions ip n are as in |(2J) but (t, z) are 
replaced by the Minkowski coordinates (t, z) adapted to 
the cavity's new rest frame, with the surface t = coin- 
ciding with 77 = 771. The nonvanishing anticommutators 
arc 



m 



{a m ,a n } = {b m ,bl\ = 5, n 



(19) 
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The Bogoliubov transformation between the Region I and 
Region III modes can then be written as 

n m 

where the coefficient matrix A = (A mn ) has the decom- 
position 

A = A^G{ m )A (21) 

and G(r/i) is the diagonal matrix whose diagonal elements 
are 

Gnnivi) = exp(ifi„77i) =: G n . (22) 

The role of G(r/i) in (|2"Tj) is to compensate for the phases 
that the modes ip m develop between rj = and r\ = rji, 
and the matrix A< = A^ 1 in (J^J) arises from matching 
Region I to Region IE at rj = r/i. Note that A is unitary 
by construction. 

Expanding A in a Maclaurin series in h as 

A = A {0) + A {1) + A (2) + 0{h 3 ) , (23) 

where the superscript again indicates the power of h, we 
obtain from (JTSJ) and ([2"T]) 

^ (0) =G(?7i), (24a) 

= G(rn)AW + (A^Gim) , (24b) 

A& = G(m)AW + (A^G( m ) + (A^G( m ) . 

(24c) 

Note that the diagonal elements of A^ arc vanishing. 
Unitarity of A implies the pcrturbative relations 

= G(r 11 )*A^ + (A^) i G(r 11 ), (25a) 

= G( m )*A^ + {A^G( m ) + (A^A^ , (25b) 
which will be useful below. 

C. Operators and vacua 

We denote the Fock vacua in Regions I and IH by | ) 
and | ) respectively. To relate the two, we mimic the 
bosonic case [TJI and make the ansatz 

|0) = Ne w | 6) , (26) 

where 

W= J2 V Pi~ a lK ( 27 ) 

p>0,g<0 

and the coefficient matrix V = (V pq ) and the normalisa- 
tion constant N are to be determined. Note that the two 
indices of V take values in disjoint sets. 



It follows from fl2a|) . fl"8|) and (J2DJ) that the creation 
and annihilation operators in Regions I and III are related 
by 

n > : a n = {j> n ,ip) = ^ a m A mn + ^ &m A mn , 

m>0 m<0 

(28a) 

n< : bl = ( ip n , ip ) = ^ a m A mn + ^ A mn . 

m>0 m<0 

(28b) 

Using (f2"6"|) and (|28ap , the condition a n | ) =0 reads 

j ]T 5 m A™« + ~ b ™ A ™n J e w | ) = . (29) 

\m>0 m<0 / 

From the anticommutators (fT9|) it follows that 

[W, a m ] = - J2 v ™ q K> ( 30a ) 

q<0 

[W,[W,a m ]] = 0. (30b) 
Using (|30]) and the Hadamard lemma, 

e' 4 Be~ A = B + [A, B] + \ [A, [A, B)} + 

(31) 

(|29|) reduces to 

A ^n V mq = - A qn (n > , q < 0) . (32) 

m>0 

A similar computation shows that the condition 6 n | ) = 
reduces to 

J2 A mn Vpm = A* pn (n < , p > 0) . (33) 
m<0 

If the block of .A where both indices are non-negative is 
invertible, Eq. (|32[1 determines V uniquely. Similarly, if 
the block of A where both indices are negative is invert- 
ible, Eq. (|33[) determines V uniquely. If both blocks are 
invertible, it can be verified using unitarity of A that the 
ensuing two expressions for V are equivalent. Working 
perturbatively in h, the invertibility assumptions hold, 
and using (f2"B"|) and ([2"4"]l we find 

V = U (1) + 0{h 2 ) (34) 

where 

V n = -4p G p = A m* G i (P > 0, 9 < 0). (35) 

We shall show in Section IIVI that the normalisation con- 
stant N has the small h expansion 

N = 1 - ^El^l 2 +0(h 3 ). (36) 
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IV. EVOLUTION OF ENTANGLED STATES 

In this section we apply the results of Section [TTT] to 
the evolution of Bell-type quantum states between the 
two cavities which are initially maximally entangled. We 
shall work perturbatively to quadratic order in h. 

Focusing first on Rob's cavity only, we write out in 
Sec. IIV Al the Region I vacuum and the Region I states 
with a single (anti-) particle in terms of Region III excita- 
tions on the Region III vacuum. In Sec. IIV Bl we address 
an entangled state where one field mode is controlled by 
Alice and one by Rob. In Sec. IIV CI we address a state of 
the type analysed in Q where the entanglement between 
Alice and Rob is in the charge of the field modes. 



A. Rob's cavity: vacuum and single-particle states 

Consider the Region I vacuum | ) in Rob's cavity. We 
shall use (|26l) to express this state in terms of Region H 
excitations over the Region H vacuum | ) . 



We expand the exponential in 
e w = l + Y^ v vA~ h \ 



as 



p,q 

+ 1 E V^afbl-afb] + 0(h 3 ). 

We denote the Region III single-particle states by 
|U) + :=4|6) 

for k > and by 

|1*>~:=&1|6) 



(37) 



(38) 



(39) 



for k < 0, so that the superscript ± indicates the sign of 
the charge. From (|37|) we obtain 



o> = | o> + $> pg \i P ) + \i q y 

~ \ E V pq V t] {l-5 pi ){l-5 qj )x 

x \i p ) + \i t ) + \i q y \i 3 y + o(h 3 ) , (40) 



where the ordering of the single-particle kets encodes the 
ordering of the fermion creation operators. It follows that 
the normalisation constant N is given by (|3"6]) . and (j2"o]l 
gives 



10) = 



i-iEi^ni^+E^i 5 *) ii«> 

P,Q P,Q 

§ E V P9 V ij( 1 - S pi)( 1 - 5 d) x 

p,q,i,j 

x |l p > + |lA + |l g )-|l.)-+0(/i 3 ). (41) 



Consider then in Rob's cavity the state with exactly 
one Region I particle, := foj, 1 ) for k < or 

|lfe) + := a' k 1 0) for k > 0. Acting on the Region I 
vacuum (|4"Tj) by (|28b[) and the Hermitian conjugate 
of (|28ap respectively, we find 



fc<0: |lfe) =Y^V m Ap k \l q ) +E^ 



rn<0 



i-iEi^i 2 ) I 1 ™) +E^(i-^ 9 )|i P ) + |i,) |i ro y 



p,q 

t \+ir \+ir 



1 E V pq V tJ (l~6 pi )(l-8 qj )(l-6 mq )(l-S mj ) |l p > + |l*> + |l ff > |li> |l m ) 



+ 0(h 3 ) , (42a) 



k > : |lfe} + = -E V P^*qk | 1p) + E A *mk 



1 - \ E i^i 2 ) 1 + E v pi c 1 - <w) 1 1 ip) + 1 1<*>~ 



+ 0{h 3 ) . (42b) 



\ E V^qV^ ( 1 - 5 pi ) (1 - ) (1 - (1 - ^i) | lr«> + | lp> + | li> + | 1q> 



p,q,i,j 
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B. Entangled two-mode states 

We wish to consider a Region I state where one field 
mode is controlled by Alice and one by Rob. Concretely, 
we take 



(43a) 



I <^it ) A R- = ^(l°fc>Al°fe)fl ± \h)l\ U )r ) . 

(43b) 

where the subscripts A and R refer to the cavity and the 
superscripts ± indicate whether the mode has positive or 
negative frequency, so that n = + for k > and k = — 
for k < 0. Furthermore, we consider the two particle ba- 
sis state of the two mode Hilbert space, corresponding to 
one excitation each in the modes k in Alice's cavity and k 
in Rob's cavity, to be ordered as in (|4*3")l . As pointed out 
in Ref. [l9| , making such a choice can lead to ambiguities 
in the entanglement. In our case, the ambiguity amounts 
to a relative phase shift of it, i.e., a sign change, in (|4"3")l . 
which does not affect the amount of entanglement. In 
other words, the states (|4"3"j) are pure, bipartite, maxi- 
mally entangled states of mode k in Alice's cavity and 
mode k in Rob's cavity. 

We form the density matrix for each of the states ((431) , 
express the density matrix in terms of Rob's Region II 
basis to order h 2 using (|4"Tj) and (|4"2j) . and take the par- 
tial trace over all of Rob's modes except the reference 
mode k. All of Rob's modes except k are thus regarded 
as environment, to which information is lost due to the 
acceleration. The relevant partial traces of Rob's matrix 
elements depend on the sign of the mode label k. For 



k > 0, corresponding to (|43a|l . we find 



Tr. fc | fc ) (0 fc | = (1 - fj~) | fc ) <0 fe | + f k |l fc ) ++ (l fc 



(44a) 



Tr. fc | 0*^1*1 = (G k + A%) | fe ) + (l fc | , (44b) 



Tr. fc |l fc ) ++ (lfe| = (1 - ft) \lkf + (l k \ + ft | 6 fc ) (6 fc | , 



(44c) 



where we have used ()25a[) and ([33]) and introduced the 
abbreviations 

-EMS/I 3 . («) 

p>0 q<0 

For k < 0, corresponding to (|43b|) . we find similarly 



Tr. fc |0 fe ) (0 fc | 
Tr^ fe | fe >-(l fe | 



(i-/ fc + )|o fe >(o fe | + /+|i fc > (l fc |, 

(46a) 

G:+A%>*)\0 k y(l k \, (46b) 

(i-/*")| i kMi*l + A~|6*)<0fc|- 

(46c) 



C. States with entanglement between opposite 
charges 



Wc finally consider the Region I state 



Xinit ) 



AR 



AO 



k)A\ L k'/R 



± II 



k' )a I l k) R 



(47) 



where the meaning of the subscripts and superscripts is as 
described for Eq. (|43|) . indicating that k > and k' < 0. 
In this state Alice and Rob each have access to both of the 
modes k and k', and the entanglement is in the charge of 
the field modes, similarly to the states considered in 

We form the reduced density matrix to order h 2 as in 
Sec. IIVB1 but now the partial tracing over Rob's modes 
excludes both mode k and mode kl . The relevant matrix 
elements take the form 

Tr^ fe , lifc'RM = / fe " 1 6 fc ) + | 0k>Y~(0 k ,\ + (b k \ 

+ (i-^-^+i^r)i^fiivr<ivr<o*i 
+ (4--i4^r)ii fe > + iM"<^r<^i 



t/<0 

+ h.c)j , (48a) 

Tr^|l fc ) ++ (lfc| = A + |6 fc ) + |6 fe ,)~(o fe ,| + (o fc | 
+ (i - ft ft + kffi I 2 ) I 1 6,' >"(0 fe< I + (i fe I 

+ (ft-\Att\ 2 )\uf\uT^\i k \ 



- (22G k G* k ,A$*A$, |6 fc > + |6 fc ,)--(l fc ,| + (l fc | 

p>0 



+ /i.e.) , 



(48b) 



+ A* kk A* k , k ,) \l k f\0 k >) (l k >\ + (0 k \ 



where in (|48c|) A kk A k , k , is kept only to order h 2 in the 
small h expansion 



AUA 



kk-^k'k' 



/m * /I (2)* 
' Ur fc' /i fc'fc' 



0(/i 3 ). 



(49) 



V. ENTANGLEMENT DEGRADATION AND 
NONLOCALITY 



We are now in a position to study the entanglement 
and the nonlocality of our states in Region IE. 
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A. Entanglement of two-mode states 



Consider the states 



btmt) ar+ and I #rit> Aii- 
in which Alice and Rob control one mode each. We shall 
quantify the entanglement by the negativity poU^ and 
the nonlocality by a possible violation of the CHSH in- 
equality mtn- 

The negativity J\f[p] is an entanglement monotone that 
quantifies how strongly the partial transpose of a density 
operator p fails to be positive. It is defined as the sum 
of the absolute values of the negative eigenvalues A of 
(1 ® T R )p, 




m = E i a i 



(50) 



A<0 



where (1 ® Tr) denotes the transpose in one of the two 
subsystems (which we have taken to be Rob without loss 
of generality). The negativity is a useful measure for 
our system because all the entangled states that it fails 
to detect are necessarily bound entangled, that is, these 
states cannot be distilled [23]] . and a system with two 
fermionic modes cannot be bound entangled. 

We work perturbatively in h. The unperturbed part 
of (1 ® Tr)pa R± has the triply degenerate eigenvalue \ 
and the nondegenerate eigenvalue — ^. In a perturba- 
tive treatment the positive eigenvalues remain positive 
and the only correction to the negativity comes from 
the perturbative correction to the negative eigenvalue. 
A straightforward computation using (|44[) and (|4l)|) shows 
that the leading correction to the negativity comes in or- 
der h 2 , and to this order the negativity formula reads 



M\Par±\ = 5 (!"/*) 



where f k := f£ + f k 



fk can be expressed as 



h= E l^- p -i| 2 l4 1)|2 



hp 



2 [Q(2k + s, 1) - Q(2k + s, E{)] h 2 



(51) 



(52) 



FIG. 2. The plot shows fk/h 2 as a function of u := 
|r7i/ln(&/a) = hTi/[45 atanh(/i/2)l , over the full period < 
u < 1. The solid curve (black) is for s — with k = ±1. 
The dashed, dash-dotted and dotted curves are respectively 
for s = t, s = | and s = |, for k = 1 (blue) above the solid 
curve and for k = — 1 (red) below the solid curve. 



with period 28 (h/ 2) atanh(/i/2), which is the proper 
time measured at the centre of Rob's cavity between send- 
ing and recapturing a light ray that is allowed to bounce 
off each wall once, fk is non-negative, and it vanishes 
only at integer multiples of the period, fk is not even in 
k for generic values of s, but it is even in fc in the limiting 
case s — in which the spectrum is symmetric between 
positive and negative charges, fk diverges at large |fc| 
proportionally to fc 2 , and the domain of validity of our 
perturbative analysis is \k\h <C 1. Plots for fc = ±1 are 
shown in Fig. [2] 

Another, potentially useful measure of entanglement 
for the states at hand would be the concurrence [25| . 
While a perturbative computation of the concurrence 
would as such be feasible, we have verified that obtain- 
ing the leading order correction would require expanding 
the states (|42[) to order h 4 . We have not pursued this 
expansion. 

We now turn to nonlocality, as quantified by the viola- 
tion of the CHSH inequality [ill, [lj] 



where 



Q(a, z) := —r Re 



a 2 ( Li 6 (z) - ^ Li 6 (z 2> 



+ Im{z)- — Li 4 (z 2 ) 
lo 



Li is the polylogarithm [24[ and 



(53) 



El :=eXP {HbJaj) =CXP U^tanh(V2) ) ' (54) 



We see from (|5ip that acceleration does degrade the 
initially maximal entanglement, and the degradation is 
determined by the function fk (|52[) . fk is periodic in t\ 



ificHSH),! < 2, 



where Schsh is the bipartite observable 



(55) 



Schsh := a • a <g> (b + b') ■ a + a' • a ® (b - b') • a . 



(56) 



a, a', b and b' are unit vectors in R 3 , and a is the vector 
of the Pauli matrices. The inequality (f55"j) is satisfied by 
all local realistic theories, but quantum mechanics allows 
the left-hand side to take values up to 2y/2. The violation 
of (|55| is hence a sufficient (although not necessary [f| 
H||) condition for the quantum state to be entangled. 

To look for violations of (|55)l . we proceed as in Q, 
noting that the maximum value of the left-hand-side in 
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the state p is given by fl2j | 

(B max ) p = 2 + fi 2 , (57) 

where pi and //2 are the two largest eigenvalues of the 
matrix U (p) — TjT p and the elements of the correlation 
matrix T = (ty) are given by £y = Tr [per j ® tjj]. In our 
scenario 

A-/* o o \ 

^(/»A/i ± )= I"/* + 0(^ 4 ), (58) 
V \-h) 

and working to order h 2 we hence find 

(Bm~) p ± R± = 2V2(l-i/ fc ) . (59) 

The acceleration thus degrades the initially maximal vi- 
olation of the CHSH inequality, and the degradation is 
again determined by the function 

B. Entanglement between opposite charges 

Wc finally turn to the entanglement between opposite 
charges in the state | xtiit ) ar S3- 

Expressing the density matrix in the Region H basis, 
tracing over Rob's unobserved modes and working pcrtur- 
batively to order h 2 , we find that the only nonvanishing 
elements of the reduced density matrix are within a 6 x 6 
block. Partially transposing Rob's subsystem replaces 
the last lines in (|48a[) and (|48b[) by their respective con- 
jugates and shifts the particle-antiparticle off-diagonals 
(|48c[) away from the diagonal. The only nonvanishing el- 
ements of the partial transpose are thus within an 8 x 8 
block, which decomposes further into two 3x3 blocks 
that correspond respectively to (|48al) and (|48b[) and the 
2x2 block 

±1 / G k G k> \A[ 1] k ,\ 2 + An,A,,\ 

2 \G^,\A^,\ 2 + At k At lk , )' 

(60) 

where the off-diagonal components are kept only to order 
h 2 in their small h expansion (|4"9")l . 

The only negative eigenvalue comes from the 2x2 
block (|60p . We find that the negativity is given by 

^] = l-iEK ) l 2 -lEI^Sr 

p/fc' p^k 

= h-\{h+fk>) + \\Et k ' -A 2 \4l\ 2 - 

(61) 

The entanglement is hence again degraded by the accel- 
eration, and the degradation has the same periodicity in 
T\ as in the cases considered above. The degradation now 
depends however on k and k' not just through the indi- 
vidual functions fk and /&< but also through the term 



proportional to | A kk , | in (|61[) : this interference term is 
nonvanishing iff k and k' have different parity, and when 
it is nonvanishing, it diminishes the degradation effect. In 
the charge-symmetric special case of s = and k = —k', 
the degradation coincides with that found in ([FT]) for the 
two-mode states (|43|) . 

VI. ONE-WAY JOURNEY 

Our analysis for the Rob trajectory that comprises Re- 
gions I, I and 11 can be generalised in a straightforward 
way to any trajectory obtained by grafting inertial and 
uniformly-accelerated segments, with arbitrary durations 
and proper accelerations. The only delicate point is that 
the phase conventions of our mode functions distinguish 
the left boundary of the cavity from the right boundary, 
and in Sec. MI Al we set up the Bogoliubov transforma- 
tion from Minkowski to Rindlcr assuming that the ac- 
celeration is to the right. It follows that the Bogoliubov 
transformation from Minkowski to leftward-accelerating 
Rindlcr is obtained from that in Sec. IIII Al by inserting 
the appropriate phase factors, A mn — > (— l)' n+n A mn , and 
in the expansions (|17[) this amounts to the replacement 
h -> —h. 

As an example, consider the Rob cavity trajectory that 
starts inertial, accelerates to the right for proper time n 
as above, coasts inertially for proper time ti and finally 
performs a braking manoeuvre that is the reverse of the 
initial acceleration, ending in an inertial state that has 
vanishing velocity with respect the initial inertial state. 
Denoting the mode functions in the final inertial state by 

ip n , and writing 

n 

we find 

Kl\ 2 - \K- n - l\ 2 \(E lE2 ) m ~ n - l\ 2 \A£l\ 2 (63) 
where E% := cxp(i7rr2/i5). For the two-mode initial states 

I $rit )ar+ and I &iit )ar- 63' the ne g ativit y and the 
maximum violation of the CHSH inequality hence read 
respectively 

^[Par±\ = |(l -Jfc), (64a) 

(ftna*)^ = 2V2"(l-# fe ), (64b) 

where 

oo 

/*= E KT 

P— — QO 

= 2 [2Q(2k + s, 1) - 2Q(2fc + 8,Ej) + Q(2k + s, E 2 ) 
- 2Q(2fc + s, EiE 2 ) + Q(2k + s, E 2 E 2 )] h 2 . 

(65) 
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v 



FIG. 3. The plot shows f k as a function of u := 
fin/ [48 atanh(/i/2)] and v := T2/(2<$) over the full period 
< u < 1 and < v < 1, for s = and fc = 1. Note 
the zeroes at u = mod 1 and at u + « = mod 1. 

The negativity in the state | xJit )aa reads 

^1 = 3 - !(/*+/*) 

+ i|£^'-l| 2 |( J B 1 E 2 ) fc - fc '-l|>«| 2 . 

(66) 

The degradation caused by acceleration is thus again 
periodic in n with period 2S(h/2) 1 atanh(/i/2), and it is 
periodic in ti with period 28. The degradation vanishes 
iff Ei = 1 or EiE% = 1, so that any degradation caused by 
the accelerated segments can be cancelled by fine-tuning 
the duration of the inertial segment, to the order h 2 in 

which we are working. A plot of f k is shown in Fig. [3] 

VII. CONCLUSIONS 

We have analysed the entanglement degradation for 
a massless Dirac field between two cavities in (1 + 1)- 
dimensional Minkowski spacctime, one cavity inertial and 
the other moving along a trajectory that consists of iner- 
tial and uniformly accelerated segments. Working in the 
approximation of small accelerations but arbitrarily long 
travel times, we found that the degradation is qualita- 
tively similar to that found in [l(| for a massless scalar 
field with Dirichlet boundary conditions. The degrada- 
tion is periodic in the durations of the individual inertial 
and accelerated segments, and we identified a travel sce- 
nario where the degradation caused by accelerated seg- 
ments can be undone by fine-tuning the duration of an 
inertial segment. The presence of charge allows however a 
wider range of initial states of interest to be analysed. As 



an example, we identified a state where the entanglement 
degradation contains a contribution due to interference 
between excitations of opposite charge. 

Compared with bosons, working in a fermionic Fock 
space led both to technical simplifications and to techni- 
cal complications. A technical simplification was that 
the relevant reduced density matrices act in a lower- 
dimensional Hilbcrt space because of the fermionic statis- 
tics, and this made it possible to quantify the entangle- 
ment not just in terms of the negativity but also in terms 
of the CHSH inequality It would further be possible 
to investigate the concurrence, although doing so would 
require pushing the perturbative low-acceleration expan- 
sion to a higher order than we have done in this paper. 

A technical complication was that when the bound- 
ary conditions at the cavity walls were chosen in an ar- 
guably natural way that preserves charge conjugation 
symmetry, the spectrum contained a zero mode. This 
zero mode could not be consistently omitted by hand, 
but we were able to regularise the zero mode by treating 
the charge-symmetric boundary conditions as a limiting 
case of charge-nonsymmetric boundary conditions. All 
our entanglement measures remained manifestly well de- 
fined when the regulator was removed. 

Another technical complication occurring for fcrmions 
is the ambiguity in the choice of the basis of the 
two-fermion Hilbert space in (|4"5)) . An alternative valid 
choice of basis is obtained by reversing the order of the 
single particle kets in (|48| , which amounts to a change of 
the signs in the off-diagonal elements of (|48a[) and (I48bl) . 
While our treatment does not remove this ambiguity, all 
of our results for the entanglement and the nonlocality of 
these states are independent of the chosen convention. 

Our analysis contained two significant limitations. 
First, while our Bogoliubov transformation technique can 
be applied to arbitrarily complicated graftings of iner- 
tial and uniformly accelerated cavity trajectory segments, 
the treatment is perturbative in the accelerations and 
hence valid only in the small acceleration limit. We 
were thus not able to address the large acceleration limit, 
in which striking qualitative differences between bosonic 
and fermionic entanglement have been found for field 
modes that arc not confined in cavities 0-0. 

Second, a massless fermion in a (1 + l)-dimensional 
cavity is unlikely to be a good model for systems realis- 
able in a laboratory. A fermion in a linearly-accelerated 
rectangular cavity in (3 + 1) dimensions can be reduced 
to the (1 + l)-dimcnsional case by separation of variables, 
but for generic field modes the transverse quantum num- 
bers then contribute to the effective (1 + l)-dimensional 
mass; further, any foreseeable experiment would presum- 
ably need to use fermions that have a positive mass al- 
ready in (3 + 1) dimensions before the reduction. It would 
be possible to analyse our (1 + l)-dimcnsional system for 
a massive fermion, and we anticipate that the mass would 
enhance the magnitude of the entanglement degradation 
as in the bosonic situation [l(J. A detailed analysis of a 
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massive fermion could become of experimental interest if 
guided by insights as to how a massive fermion might be 
confined to a cavity in a concrete laboratory setting. 

We started this paper by emphasising that a cavity lo- 
calises the quantum degrees of freedom in the worldtubc 
of the cavity, and our assumption of inertial initial and 
final trajectory segments localises the acceleration effects 
in a finite interval of the cavity's proper time. We should 
perhaps end by emphasising that we are not attempt- 
ing to localise measurements of the field at more precise 
spacetime locations within the worldtubc of the cavity, 
and we are hence not proposing cavities as a fundamental 
solution to the open conceptual issues of a quantum mea- 
surement theory in relativistic spacetime [271 ] . A cavity 
can however reduce the measurement ambiguities from, 



say, megaparsecs to centimetres, which may well suffice 
to resolve the conceptual issues in specific experimental 
settings of interest, gedanken or otherwise. 
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